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Rate  Constant  vs.  Time  Near  Cloud  Base 


RADIATION  DOSES  FROM  FLYING 
THROUGH  NUCLEAR  DEBRIS  CLOUDS 


INTRODUCTION 

Estimates  of  radiation  doses  to  aircrews  flying  through  nuclear  debris 
clouds  are  of  critical  interest  to  USAF  mission  planners.  The  probability 
of  such  an  event  occurring  is  high  following  a  nuclear  attack.  Since  real  data 
are  lacking,  prediction  must  depend  upon  mathematical  modeling  of  the  cloud 
environment  and  fly through. 

Taboada  et  al.  (1)  recently  developed  a  computer  model  which  predicts  doses 
from  flythroughs.  The  model,  which  will  be  called  CASSIE  in  this  report,  is 
based  upon  the  dust  environment  model  CASSANDRA  (2),  which  in  turn  uses  the 
benchmark  DELFIC  (3)  code  as  the  basis  of  its  cloud  rise  and  fallout  dynamics. 
The  CASSIE  model  is  user  friendly;  it  supplies  default  values  and  descriptions 
of  parameters  when  requested  and  has  all  the  input  flexibilities  of  the  DELFIC 
and  CASSANDRA  codes. 

The  CASSIE  model  uses  the  CASSANDRA  code  to  calculate  the  dust  density  at 
various  points  along  the  trajectory  of  an  aircraft  flying  through  a  nuclear 
debris  cloud.  A  gamma- radiation  dose  rate  is  calculated  from  the  dust  densities 
and  integrated  over  time  to  obtain  the  total  radiation  dose  that  an  aircrew 
would  experience. 

The  total  radiation  dose  is  assumed  to  originate  from  2  sources:  (1)  the 
immersion  radiation  dose,  which  is  a  result  of  the  aircraft  being  immersed  In  a 
radioactive  cloud,  and  (2)  the  onboard  radiation  dose,  which  results  from  radia¬ 
tion  exposure  of  the  aircrews  to  dust  particles  trapped  inside  the  aircraft 
cabin  and  filtering  system.  The  immersion  rate  exists  only  while  the  aircraft 
is  inside  the  radioactive  debris  cloud,  whereas  the  onboard  radiation  is  a  con¬ 
tinuing  hazard. 

The  purpose  of  this  report  is  to  use  the  CASSIE  model  to  calculate  the 
radiation  doses  to  aircrewmembers  for  various  conditions  of  flythrough.  Due  to 
the  long  program  running  time,  the  model  would  have  little  use  in  an  operational 
setting.  The  intent  of  this  report  is  to  supply  results  of  the  model  in  a  form 
which  will  allow  planners  to  estimate  in  advance  realistic  and  worst  case  doses 
for  flythrough  scenarios. 


FLYTHROUGH  MODEL 

CASSIE  provides  for  changing  flight  and  aircraft  parameters  easily.  A 
scenario  is  needed  which  will  allow  an  analysis  of  radiation  dose  dependence  on 
height  of  flythrough,  time  of  flight,  and  detonation  yield.  For  the  purposes  of 
this  report,  the  following  flight  scenario  was  used. 

An  aircraft  flies  from  a  point  40  km  (24  mi)  from  ground  zero  on  one  side 
of  a  detonation  to  40  km  on  the  opposite  side,  passing  directly  over  ground 
zero.  The  flight  duration  i3  500  s,  which  corresponds  to  a  velocity  of  311 
knots.  A  height  H  is  maintained  throughout  the  flight  and  the  aircraft  passes 


over  ground  zero  at  t  seconds  after  detonation.  This  path  is  expressed  in  the 
form  (-40  km,  0,  H,  t-250)  -  (40  km,  0,  H,  t+250).' 


The  aircraft  parameters  are  chosen  to  resemble  those  of  a  KC-135  aircraft. 
The  model  approximates  the  cabin  as  a  cylinder,  chosen  to  be  32.05  m  (107  ft)  in 
length  and  2.74  m  (9  ft)  in  radius.  The  filtration  system  is  arbitrarily  chosen 
to  allow  50$  of  the  dust  particles  taken  onboard  to  enter  the  cabin,  with  the 
remainder  staying  on  the  filter.  Once  aboard,  the  dust  particles  do  not  exit. 
The  distance  from  the  crewmembers  to  the  filter  is  arbitrarily  chosen  to  be  250 
cm  (100  in.).  The  air-mass  flow  into  the  aircraft  is  67.5  kg/min  (150  lb/min). 
The  dose  received  by  aircrewmembers  will  be  dependent  on  the  choice  of  these 
parameters,  which  are  discussed  in  Appendix  A.  If  cloud  penetration  is  unavoid¬ 
able,  Figure  3  shows  that  it  is  most  advantageous  to  fly  through  as  high  in  the 
cloud  as  possible.  This  will  allow  the  aircraft  to  be  exposed  to  the  lowest 
possible  dose  within  the  cloud.  A  flythrough  beneath  the  cloud  may  result  in  a 
lower  dose,  but  there  is  the  danger  of  the  aircraft  encountering  large  parti¬ 
cles. 


CLOUD  DIMENSIONS 

The  size  of  the  nuclear  debris  cloud  depends  on  a  number  of  factors.  The 
parameters  used  in  the  analysis  include: 


(1) 

A  surface 

burst,  with  50$ 

of  the  available  energy  expended  in  fission. 

(2) 

Ground  zero  at  670  m  (2211 

ft)  above  mean  sea  level. 

(3) 

No  wind. 

(4) 

Soil  type 

is  siliceous. 

The  debris  cloud  is  modeled  throughout  its  creation,  expansion,  and  cool¬ 
ing.  At  some  point  the  model  considers  the  turbulence,  upward  rise,  and  expan¬ 
sion  to  be  negligible  and  stabilizes  the  cloud  dimensions  at  current  values. 
These  calculations  are  reported  in  Table  1  and  the  cloud  dimensions  are  shown 
in  Figures  1  and  2.  Stabilization  is  relative;  in  actuality,  the  cloud  will 
continue  to  grow  slowly  as  a  result  of  buoyancy,  turbulence,  and  wind  patterns. 


TABLE  1.  CLOUD  DIMENSION  DATA  AT  STABILIZATION 


Yield  (MT) 

Stab . 
Time(s) 

Cloud 

Base(m) 

Cloud 

Top 

Cloud 

Radius 

0.003 

312 

3700 

5200 

1100 

0.01 

423 

5300 

7400 

1800 

0.03 

661 

6700 

10500 

2900 

0.1 

611 

8400 

12400 

5400 

0.3 

720 

9500 

14900 

8500 

1  . 
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Figure  1.  Base  and  top  heights  of  clouds  vs.  yield 


These  calculations  are  based  on  the  DELFIC  and  CASSANDRA  portions  of  the 
model.  The  cloud  top  and  bottom  heights  agree  well  with  values  reported  by 
Glasstone  and  Dolan  (4).  Cloud  radii  appear  to  be  larger  by  as  much  as  a  factor 
of  1.4  than  those  reported.  Part  of  the  disagreement  may  be  explained  by  the 
ambiguous  nature  of  the  cloud  dimensions.  Part  may  be  explained  by  the  fact 
that  DELFIC  was  intended  as  a  ground  fallout  code,  and  ground  fallout  is  rela¬ 
tively  insensitive  to  cloud  radius  or  base  height  choices  (7).  This  is  a  draw¬ 
back  since  flythrough  dose  is  approximately  inversely  proportional  to  cloud 
radius  in  this  model. 


DUST  CLOUD  CHARACTERISTICS 

In  the  model,  the  mass  of  soil  entrained  by  the  burst  is  dependent  on  the 
height  of  burst  above  ground,  type  of  soil,  and  explosive  yield.  A  log-normal 
distribution  of  particle  sizes  is  assumed,  and  the  dust  is  partitioned  into  100 
size  classes.  At  initial  time  the  soil  burden  is  distributed  uniformly  through¬ 
out  1  the  cloud.  As  the  cloud  grows  with  time,  the  soil  debris  will  be  lofted 
with  the  cloud  at  a  rate  dependent  on  particle  size.  The  gamma- radiation  rate 
at  each  point  in  the  cloud  is  proportional  to  the  dust  mass  concentration  at 
that  point.  The  latter  assumption  is  considered  in  detail  in  Appendix  3. 

One  way  of  looking  at  the  vertical  distribution  of  dust  is  to  conduct  fly- 
throughs  at  various  heights.  Figure  3  shows  the  dependence  of  total  dose  accu¬ 
mulated  during  the  flythrough  scenario  as  height  is  varied  for  3  different 
times.  For  T-500  s  the  flythrough  starts  250  s  and  ends  750  s  after  burst.  The 
aircraft  passes  through  the  entire  cloud  (radius  -  13.8  km  (8.4  mi))  before 
stabilization  is  reached  at  783  3  (Table  1).  The  dose  is  nearly  constant  for 
flythroughs  passing  through  the  cloud.  As  the  aircraft  flies  at  heights  lower 
than  the  cloud  base,  the  dose  received  is  relatively  much  lower;  only  large 
particles  have  fallen  out  to  those  levels.  At  later  times,  the  dust  activity 
has  decayed  as  well  as  fallen.  The  dose  received  during  flights  below  the  cloud 
base  is  relatively  more  substantial. 

The  flythrough  at  the  cloud  base  receives  the  largest  dose.  In  the  follow¬ 
ing  analyses,  this  will  be  considered  as  the  worst  case  scenario. 


DOSE  AT  THE  END  OF  FLYTHROUGH 

The  dose  from  a  flythrough  is  strongly  dependent  on  the  time  of  flythrough 
and  the  yield  of  the  detonation.  Figure  4  shows  the  dose  accumulated  by  the  end 
of  the  flythrough  at  the  cloud  base  vs.  yield  at  1000  s.  Although  the  data  on 
this  graph  (Fig.  4)  do  not  all  lie  on  a  straight  line,  the  data  with  yields 
greater  than  0.1  MT  strongly  suggest  a  power  relationship. 

Figure  5  3hows  the  dependence  of  dose  at  the  end  of  flythrough  on  time  of 
flythrough.  The  time  plotted  is  the  time  after  detonation  that  the  aircraft 
passes  over  ground  zero.  Regression  lines  to  these  data  are  also  shown.  The 
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Figure  5.  Dose  vs.  time  for  a  10  MT  yield. 
Flythrough  (-40km,  0,  H,  t-250)-(«»0km,  0,  H,  t+250) 


regression  fit  to  the  flythrough  at  an  altitude  of  2  km  (1.2  mi)  does  not  in¬ 
clude  the  data  for  the  four  earliest  times.  The  flythrough  at  19-km  (11.6  mi) 
altitude  represents  a  flight  through  the  base  of  the  cloud  while  the  flythrough 
at  28-km  (17.1  mi)  altitude  is  close  to  the  top  of  the  cloud.  The  power  depend¬ 
ence  of  dose  on  time  is  greater  than  the  t"1 *2  (Way-Wigner  (5))  dependence  built 
into  the  model.  This  is  a  result  of  the  dust  "falling  out"  to  lower  levels. 


Figures  6  and  7  show  the  dependence  of  dose  on  time  for  flythroughs  near 
the  cloud  top  and  bottom,  respectively.  These  graphs  show  that  the  exponent  in 
the  time  dependence  is  nearly  independent  of  yield.  A  multiple  regression  on 
data  for  times  between  .5  h  and  3  days,  for  yields  above  .1  MT,  and  for  fly¬ 
throughs  near  the  cloud  base  produces  the  approximate  relationship 

Dose  -  2.34  W1*®  t  -1-53  ( 1 ) 

where  W  is  in  megatons  and  t  is  in  hours.  This  time  dependence  compares  favor¬ 
ably  with  the  t“  dependency  observed  in  actual  flythroughs  (6). 


POST  FLYTHROUGH  DOSE 

The  radiation  from  onboard  sources  is  insignificant  (3S  of  immersion  dose) 
while  the  aircraft  is  immersed  in  the  cloud.  However,  since  the  onboard  dust 
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Figure  6.  Dose  vs.  time  near  cloud  top. 
Flythrough  (-40  km,  0,  Ht.  t-250)-(40  km,  0,  Ht,  t+250) 
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Figure  7.  Dose  vs.  time  near  cloud  base. 
Flythrough  (-40  km,  0,  Hb,  t-250)-(-40  km,  0,  Hb,  t+250) 


particles  are  carried  with  the  aircraft,  the  radiation  from  this  source  can 
become  significant  for  extended  missions.  The  gamma-radiation  rate  decays  as: 

R(t)  -  R1  t'1*2  (2) 

where  R(t)  is  the  radiation  dose  rate  in  rads/hour  at  time  t  hours  after  detona¬ 
tion  and  R^  is  the  dose  rate  constant,  equal  to  the  dose  rate  at  one  hour  after 
detonation  (5).  The  evaluation  of  R1  is  based  on  the  actual  values  of  the  total 
onboard  dose  rates  at  the  corresponding  cloud  exit  times  of  the  aircraft.  The 
calculated  dose  rate  constants  for  trajectories  near  the  cloud  base  and  near  the 
cloud  top  are  given  in  Figures  8  and  9.  A  multiple  regression  fit  to  the  data 
shows  that  the  dose  rate  constant  can  be  approximated  by 

R1  -  1  .44  W55  t  "*33  (3) 

where  W  is  yield  (megatons)  and  t  is  time  (hours)  after  detonation  that  the 
plane  reaches  ground  zero. 

The  total  onboard  radiation  dose  can  be  obtained  by  integrating  equation 
(2)  over  the  total  time  of  the  mission: 


D(t)  -  jf  R1  t'1*2  dt  -  5R1  (ta  “°-2  -  t  ~0'2) 
a 


(4) 


where  D(t)  is  the  dose  in  rads  for  the  mission  at  time  t  hours  after  detona¬ 
tion  and  ta  is  the  time  after  detonation  that  the  aircraft  exits  the  debris 
cloud. 


Figure  8.  Rate  constant  vs.  time  near  cloud  top. 
Flythrough  (-40  km,  0,  Ht,  t-250)-(40  km,  0,  Ht,  t+250). 
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Figure  9.  Rate  constant  vs.  time  near  cloud  base. 
Flythrough  (-40  km,  Hb,  t-250)-(40  km,  0,  Hb,  t+250). 


CONCLUSIONS 

The  problem  of  radiation  exposure  to  aircrews  during  or  Immediately  after  a 
nuclear  attack  is  critical  for  USAF  mission  planners.  The  best  way  to  reduce 
radiation  exposure  is  for  the  aircraft  to  avoid  the  nuclear  debris  clouds.  How¬ 
ever,  the  size  of  the  cloud  or  the  number  of  clouds  may  preclude  this  option. 

The  CASSIE  model  will  calculate  the  radiation  doses  to  aircrews  flying 
through  the  debris  clouds  from  a  few  seconds  after  detonation  to  time  of  cloud 
disbursement.  However,  the  residual  nuclear  radiation  will  not  be  the  major 
concern  for  early  times.  The  initial  radiation,  thermal  effects,  and  the  prob¬ 
lems  of  flying  through  clouds  of  large  debris  particle  sizes  will  be  the  major 
considerations.  Thus,  the  results  of  this  report  should  not  be  considered  as 
feasible  for  flythrough  times  before  cloud  stabilization. 

Radiation  doses  of  aircrews  flying  through  the  nuclear  debris  clouds  under 
various  conditions  have  been  calculated  by  the  use  of  the  CASSIE  model  for  an 
aircraft  moving  with  a  constant  speed  of  about  311  knots.  The  radiation  dose 
for  an  aircraft  moving  at  a  different  speed  will  change  in  an  essentially  in¬ 
verse  relationship  with  the  speed  of  the  aircraft.  If  an  aircraft  penetrates 
the  cloud  at  a  slower  speed,  the  aircraft  will  be  exposed  to  the  radiation  for  a 
longer  time  and  thus,  the  dose  will  increase.  The  adjusted  dose,  D1 ,  is  calcu¬ 
lated  by  the  equation: 


D,  -  (311/v)  D 


(5) 


where  v  is  the  actual  speed  of  the  aircraft  in  knots  and  D  is  the  predicted 
dose  that  the  aircraft  would  have  had  if  it  traveled  at  a  speed  of  311  knots. 
The  error  introduced  by  this  approximation  is  less  than  6%  for  a  range  of  speeds 
of  200  to  700  knots  at  a  flythrough  time  of  1  h.  For  later  times  or  for  speeds 
closer  to  311  knots,  this  approximation  is  even  better. 

To  illustrate  the  application  of  the  results  of  this  study,  consider  the 
following  problem.  An  estimate  of  the  radiation  dose  is  desired  for  an  aircraft 
flying  through  the  middle  of  a  debris  cloud  at  the  cloud  base  from  a  0.5-MT 
surface  burst  with  a  flythrough  speed  of  250  knots  at  a  time  of  2000  s  after 
detonation  enroute  to  a  target  that  requires  12  h  for  mission  completion.  The 
steps  for  this  calculation  are  as  follows: 

(1)  For  the  moment,  consider  the  aircraft  as  having  a  speed  of  311 
knots.  Calculate  the  radiation  dose  while  immersed  in  the  debris  cloud  from 
equation  (1): 

D  -  (2. 31*)  (.5)  0,48  (2000/3600)“’ *53  (6) 
D  «  4.1  rads 

(2)  Find  the  onboard  radiation  dose  rate  constant  from  equation  (3): 

R1  -  (1.44)  (.5)‘55  (2000/3600)"° *33  (7) 
-  1.2  rads/hour 


(3)  Find  the  onboard  dose  for  a  12-h  mission  from  equation  (4): 

D( 12)  -  (5)  (1.2)  [(2000/3600)’* 2  -  (12)“*2]  (8) 

-  3.1  rads 

At  this  point  multiply  by  the  aircraft  correction  factor  if  known;  see  Appendix 
A. 

(4)  Calculate  the  total  dose  that  the  aircrew  would  experience  by  adding 
the  results  of  steps  (1)  and  (3): 

Total  dose  -  7.2  rads 

(5)  Adjust  the  dose  to  account  for  the  speed  of  the  aircraft  by  using 
equation  (5) : 

D1  -  (31  1/250K7.2)  -  9.0  rads 

Thus,  the  predicted  gamma -radiation  dose  to  the  aircrew  is  9.0  rads  for  the 
mission.  For  50  runs  in  the  range  .1  MTOKIO  MT  and  30  min<T<  1  day,  the  stan¬ 
dard  deviation  of  these  predictions  from  the  CASSIE  predictions  is  16%.  Two 
runs  deviated  more  than  30%  from  the  CASSIE  results. 


If  the  aircraft  flies  through  more  than  one  cloud,  the  radiation  dose  of 
each  cloud  oust  be  calculated  separately.  The  total  dose  is  the  sum  of  all  the 
individual  doses. 
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CORRECTION  FOR  AIRCRAFT  PARAMETERS 


To  obtain  the  dose  prediction  outlined  in  the  conclusions,  the  following 
aircraft  parameters  were  arbitrarily  chosen: 

CL(cabin  length)  -  3205  cm  (1300  in.) 

CR(cabin  radius)  -  274  cm  (110  in.) 

DF(distance  of  crewmember  to  filter)  -  250  cm  (100  in.) 

FR(air-flow  rate  into  cabin)  -  67.5  kg  (150  lb)/min 
PF(fraction  of  dust  not  trapped  by  filter)  -  .5 

Given  actual  measurements  of  these  5  parameters,  a  correction  factor  for  a 
particular  type  of  aircraft  may  be  calculated  and  the  carryalong  dose  modified 
as  follows: 

1.  Calculate  the  correction  factor  for  cabin  dimensions,  CDM: 

CDM  -  147060[CL-L0G(4CR2+a2)  ♦  4CR*ATN(CL/2CR)  -  2CL'L0G(CL)]/(CL*CR2) 

2.  Calculate  the  correction  factor  for  crew  distance  from  the  filter,  CF: 

CF  -  62500/D2 

3.  Calculate  the  correction  factor  for  air  flow  into  the  cabin,  CFR: 

Cj.R  -  FR/150 

4.  Calculate  the  correction  factors  for  the  fraction  of  dust  allowed  to  pass 
into  the  cabin  through  the  filter,  CpF1  and  0pF2: 

CpFi  -  2PF 
CpF2  -  2(1-PF) 

5.  Calculate  the  combined  correction  factor,  C: 

C  -  0.3'CpF1  CDM  CFR  +  0 .7 'CpF2  CFR  CF 

The  carryalong  dose  (result  of  step  3  in  the  conclusions)  is  multiplied  by 
this  combined  correction  factor  to  obtain  the  corrected  carryalong  dose.  The 
correction  factor  is  based  on:  (a)  analytical  considerations  discussed  by 

Taboada  et  al.(l),  and  (b)  the  empirical  observation  that,  given  the  arbitrary 
parameters  originally  chosen,  the  CASSIE  model  calculates  that  30. %  of  the 
carryalong  dose  is  due  to  dust  trapped  in  the  cabin  while  70. %  is  due  to  dust 
trapped  in  the  filter.  The  correction  factor  thus  contributes  no  deviation  of 
hand  calculated  results  from  CASSIE  results. 


APPENDIX  B 


EVALUATION  OF  MODEL  ASSUMPTIONS 

In  the  CASSIE  model,  the  dose  rate  for  a  flythrough  is  calculated  by  the 
following  steps. 

(1)  Divide  the  flight  path  into  a  sequence  of  target  points  (x1 ,  y1 ,  z1 ,  t1 ) , 
(x2,  y2,  z2,  t2),...(xn,  yn,  zn,  tn),  each  4  km  apart. 

(2)  Calculate  the  dust  concentration  at  each  target  point  for  each  of  100  dust 
sizes  using  the  CASSANDRA  code. 

(3)  Calculate  the  dose  rate  at  each  target  point  assuming: 

(a)  The  rate  contributed  by  a  volume  of  space  is  proportional  to  the  mass 
of  dust  contained  therein. 

(b)  When  the  dose  rate  is  calculated  at  a  particular  target  point,  the 
dust  concentration  is  uniform  and  equal  to  the  concentration  at  the  target 
point. 

(c)  The  cloud  is  considered  to  be  infinite  in  extent  for  the  purpose  of 
integrating  the  dose  rate  for  a  target  point. 

(M)  Approximate  the  dose  integration  over  the  flight  path  by  assuming: 

(a)  Constant  velocity 

(b)  The  dose  rate  varies  linearly  between  target  points. 

Each  of  these  operations  contributes  some  error.  An  analysis  of  the  extent 
of  some  of  them  follows. 


Frequency  of  Target  Points 

A  variety  of  flight  paths  were  evaluated  by  taking  points  much  closer 
together  than  4  km  (2.4  mi).  The  differences  between  these  doses  and  the 
original  dose  were  less  than  5%.  Thus,  it  was  concluded  that  the  4  km  (2.4  mi) 
interval  for  target  points  was  a  good  compromise  between  accuracy  and  computer 
time. 


Integration  to  get  Gamma  Rate 


Taboada  et  al.  (1)  calculated  the  immersion  dose  rate  D  at  a  point  (r^) 
and  time  (t)  by  the  following  integral: 


Dfrp 


100 

t"1'2  K,jK(J)/  p(r,t,j)  exp  (-K3)]r  -  r,  I )  dV 


(B-1) 


The  space  dependent  portion  of  the  gamma  rate  at  any  point  for  each  dust 
size  category  is 


D  (rifj)  -  Constant’  /  p(r  ,j)  exp  (-K,  lr  -  r^ I ) 


(B-2) 


Mir  lr  -  r. 


where  is  the  i  th  target  point  location,  p(r)  is  the  dust  density  at  loca¬ 
tion  r,  -  6.767  x  10*"3  meters  -1  is  the  atmospheric  absorption  of  a  1- 
MeV  gamma  photon,  and  j  is  the  dust  size  category.  The  integral  is  taken 
over  the  volume  of  the  cloud.  If  a  coordinate  transformation  is  made  to  put 
the  origin  at  the  target  point  and  spherical  coordinates  are  used,  this 
equation  becomes 


D  (r.,,j)  -  Constant  ///  p(s,J)  exp  (-K,  s)  sine  d0d$ds 

TT  3 


(B-3) 


If  the  assumptions  in  steps  3(b)  and  3(0)  are  made,  then  p(s)  -  pCrp.  Then, 
p(s)  is  constant  for  purposes  of  evaluating  the  integral.  Also,  the  cloud  is 
spherically  symmetric  and  infinite  in  size.  Thus,  the  integral  becomes 


D(ritj)  -  Constant  p(rifj)  /  exp  (-K^  s)  ds 

o  J 


D(r^  ,J  )  -  Constant  p(ritj)/Kg 


(B-4) 


(B-5) 


Two  objections  may  occur  at  this  point.  First,  the  dust  density  does  not 
continue  indefinitely.  This  is  not  a  serious  problem  because  only  1$  of  the 
dose  rate  calculated  by  this  approximation  is  due  to  any  activity  beyond  a  680 
m  (2244  ft)  sphere  about  r1 .  That  is, 


680  °> 

£  expC-K^  s)ds  »  .99  fgexp(~K^  s)ds 


(B-6) 


A  second  objection  is  that  p(r)  is  not  constant  enough  even  within  a  680  m 
(2244  ft)  sphere  to  make  a  good  approximation.  To  estimate  the  error  incurred 
from  the  constant  density  approximation,  an  alternative  integration  scheme  was 
performed,  using  the  following  symmetry  argument.  Consider  the  target  point  (x1, 
yi*  zi»  ti^’  po11  the  du3t  concentrations  at  7  points:  (x^,  y^,  zi(  t,),  (xt  ± 
200m,  yif  z^,  tj),  (x^,  yj  ±  200m,  zif  tj)  and  (xi(  y,,  ±  200m,  tj^).  Assume 
the  dust  concentration  is  uniform  within  a  100  ra  (327  ft)  sphere  around  the 
target  point.  This  region  would  encompass  one-half  the  dose  rate  in  the  previ¬ 
ous  scheme.  Divide  the  remaining  space  into  6  "congruent"  (infinite)  regions, 
each  with  axis  of  symmetry  about  the  line  through  (x^,  y^,  z and  the  respec¬ 
tive  polling  point.  This  can  be  pictured  by  imagining  the  8  points  of  a  cube, 
(xj  ±  d,  y^  ±  d.Zj^  ±  d),  on  a  hollow  cantaloupe  and  slicing  between  them.  Of 
course,  this  cantaloupe  is  infinite  in  radius.  The  intent  of  this  division  is 
to  Justify  a  weighting  of  the  polled  points.  The  dose  rate  due  to  one  of  the 
outside  regions  is  one-twelfth  of  the  rate  calculated  over  all  space  if  the  dust 
concentration  were  uniformly  equal  to  that  of  the  polled  point.  In  like  manner, 
the  interior  sphere  receives  a  one-half  weighting. 


'•V • 'V'-V'V 


The  sequence  of  polling  points  (x1 ,  y1 ,  z.  ♦  200,  tj) .  (xn,  yn,  zn  ♦ 

200,  tn)  creates  a  flight  path  200  m  (660  ft;  above  the  original.  The  dose 
received  from  the  upper  outside  region  can  be  evaluated  by  running  the  original 
program  and  taking  one-twelfth  of  the  resulting  dose.  The  same  is  true  of  (xif 
yA  ±  200,  Zj,  tt)  and  (xA ,  yA ,  -  200,  t^).  The  other  2  outside  integration 

regions  (Xj^  ±  200,  yi#  z. ,  t^)  form  flight  paths  along  the  original  path  but 
shifted  in  time  by  200  m  (660  ft)/velocity ,  which  is  a  negligible  amount.  It  is 
assumed  that  the  starting  and  ending  points  are  outside  the  cloud. 

Thus,  the  dose  from  this  integration  scheme  is  a  linear  function  of  5  runs 
of  the  original  program: 

New  dose  estimate 
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The  doses  predicted  by  CASSIE  along  the  (x,,  y.  ±  200,  z1 ,  t.,)  - 
(xn,  yn  ±  200,  zn,  tn)  flight  paths  do  not  differ  from  the  prediction  along  the 
original  path.  As  can  be  seen  from  Figure  3,  doses  along  the  (x1 ,  y-j.z.  ±  200, 
t1^xn»  ynf  zn  ±  200»  tn^  Patha  may  differ’  by  25$  from  the  prediction  along  the 
original  path  but  both  will  not  differ  in  the  same  direction  from  the  original. 
The  choice  of  one  flight  path,  to  differ  by  25$  while  the  other  agrees  with  the 
original  is  a  worst  case  and  results  in  25/12  or  about  2$  error.  We  take  this 
as  evidence  that  the  original  integration  scheme  is  sufficiently  accurate. 


Calculation  of  Dust  Densities 

Results  of  CASSANDRA  agree  well  with  dust  concentrations  measured  in  the 
Dial  Pack  tests  using  conventional  explosives.  There  are  no  measurements  of 
airborne  dust  from  detonations  in  the  megaton  range  with  which  to  compare  the 
model’s  results.  The  cloud  size  discrepancy  in  the  megaton  range  has  already 
been  noted.  The  lack  of  experimental  data  and  other  models  for  comparison  makes 
the  sources  and  amounts  of  error  difficult  to  assess. 

Atmospheric  Conditions 

The  atmospheric  winds  are  assumed  to  be  zero.  If  nonzero  winds  are  di¬ 
rected  horizontally,  the  cloud  will  be  translated  as  a  unit  if  the  winds  are 
uniform  and  distorted  if  a  wind  shear  exists.  A  wind  shear  along  the  direction 
of  the  flight  path  will  cause  little  difference  in  dose  while  a  perpendicular 
wind  shear  will  produce  a  lower  dose. 

The  effects  of  precipitation  and  large  vertical  air  currents  as  in  a 
thunderstorm  are  not  considered.  Precipitation  will  cause  the  dust  particles  to 
fall  out  earlier  and  thus,  reduce  doses  from  flythroughs. 
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